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Abstract. Perfect fluid tori with uniform distribution of the specific angular momentum orbiting the Kerr-de Sitter black holes 
or naked singularities are studied. Closed equipotential surfaces corresponding to stationary toroidal discs are allowed only 
in the spacetimes admitting stable circular geodesies. The last closed surface crosses itself in the cusp(s) enabling outflow(s) 
of matter from the torus due to the violation of hydrostatic equilibrium. The repulsive cosmological constant, A > 0, implies 
the existence of the outer cusp (with a stabilizing effect on the tori because of excretion, i.e., outflow of matter from the torus 
into the outer space) and the strong collimation of open equipotential surfaces along the rotational axis. Both the effects take 
place nearby so-called static radius where the gravitational attraction is just balanced by the cosmic repulsion. The plus-family 
discs (which are always corotating in the black-hole backgrounds but can be counterrotating, even with negative energy of the 
fluid elements, in some naked-singularity backgrounds) are thicker and more extended than the minus-family ones (which are 
always counterrotating in all backgrounds). If parameters of the naked-singularity spacetimes are very close to the parameters 
of extreme black-hole spacetimes, the family of possible disc-like configurations includes members with two isolated discs 
where the inner one is always a counterrotating accretion disc. Mass estimates for tori with nonrelativistic adiabatic equation 
of state give limits on their central mass-density, for which the approximation of test fluid is adequate. 
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1. INTRODUCTION 

Disc-like structures orbiting black holes seem to play an important role in a wide range of astrophysical phenomena 
including the most energetic processes in the Universe connected with quasars, or the formation of probably the largest 
disc structures at all-galactic discs. On the other hand, various cosmological observations indicate convincingly that 
in the framework of inflationary cosmology a non-zero, although very small, vacuum energy density, i.e., a relic 
repulsive cosmological constant, A > 0, or some similarly acting new kind of fields called quintessence, has to be 
invoked in order to explain the dynamics of the recent Universe [1, 2]. Both possibilities are often referred as a 
dark energy in the Universe. It is well known that the repulsive cosmological constant strongly influences expansion 
of the Universe, leading finally to an exponentially accelerated stage. The paper shows, surprisingly enough, that 
the repulsive cosmological constant could be relevant also in astrophysical processes like the disc accretion onto a 
supermassive black hole and the formation of the largest disc structures in the Universe. 

Recent data from a wide variety of independent cosmological tests give the current value of the vacuum energy 
density to be [3] 

Pvac(0) ~ 0- 7 3Pcrit(0); (1) 

where the present value of the critical energy density p cr i t (o) is related with the Hubble parameter Ho by 1 

PcdtfO) = // = 100/ ! kms- 1 Mpc- 1 . (2) 

Taking value of the dimensionless parameter /i«0.7, we obtain for the relic repulsive cosmological constant its current 
value 

A = 87rp vac(0) w 1.3 x 1(T 56 cm" 2 . (3) 



1 Geometric units c = G = 1 are used hereafter. 
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Basic properties of geometrically thin (accretion) discs with low accretion rates and negligible pressure are given 
by the circular geodesic motion in the black-hole backgrounds [4], while for geometrically thick discs with high 
accretion rates and pressure being relevant they are determined by equipotential surfaces of test perfect fluid rotating 
in the background (see, e.g., [5]). The equipotential surfaces can be closed or open. Moreover, there is a special class 
of critical surfaces self-crossing in a cusp(s), which can be either marginally closed or open. The closed — toroidal — 
equipotential surfaces determine stationary configurations (tori). The fluid can fill any toroidal equipotential surface; 
at the surface of the torus pressure vanishes, but its gradient is non-zero [15]. On the other hand, the open equipotential 
surfaces are important in dynamical situations, e.g., in modeling of jets [20, 21], as the matter can flow along the open 
surfaces. The critical — marginally closed — equipotential surfaces W cr i t are important in the theory of thick accretion 
discs, because accretion onto the black hole through a cusp of the equipotential surface, located in the equatorial 
plane, is possible due to a small overcoming of the critical surface by the surface of the disc (Paczyhski mechanism). 
Accretion is thus driven by a violation of the hydrostatic equilibrium, rather than by viscosity of the accreting matter 
[15]. 

The equatorial circular motion of test particles (Keplerian motion) and its relevance for thin accretion discs were 
studied thoroughly in the case of Schwarzschild-de Sitter (SdS) geometry in [6], and in the case of Kerr-de Sitter 
(KdS) geometry in [7], even if some important remarks can also be found in the work of Carter or Demianski [8, 9]. It 
was shown that thin discs have, besides the standard inner edge at the inner marginally stable circular orbit, an outer 
edge at the outer marginally stable circular orbit located slightly under the static radius of the given spacetime. The 
static radius is an unstable circular orbit where the gravitational attraction of the black hole is just compensated by the 
cosmic repulsion and where the static geodesic observer, i.e., the geodesic observer with only the time-component of its 
4- velocity being non-zero, resides. Moreover the accretion efficiency, 77 =£ ms ( ) — £ ms (i), given by the difference of the 
specific energies at the outer and the inner marginally stable orbits, is smaller when compared with the asymptotically 
flat geometries, especially due to the fact that E ms ( j < 1. 

Equilibrium toroidal configurations of barotropic perfect fluid orbiting in the Schwarzschild-de Sitter background 
were studied in [10]. An analysis of the structure of equipotential surfaces revealed the existence of another (outer) 
critical point (the so-called cusp) in which one of the equipotential surfaces is self-crossing. As the discussed 
equilibrium of fluid tori is the hydrostatic equilibrium, the outer cusp, when it is located on the marginally closed 
equipotential surface, determines the outer edge of the disc, exactly in the same way as the inner cusp determines 
the inner edge [15]. Moreover, the outflows through the outer cusp can stabilize the accretion discs against so-called 
runaway instability [23], as shown by Rezzolla et al. [11]. When the critical surface with the inner cusp is open or 
does not exist, while the critical surface with the outer cusp is marginally closed, such a configuration corresponds to 
a completely new type of discs called excretion disc. More detailed description of possible equilibrium configurations 
of barotropic fluid in the backgrounds with A > is given in the Sec. 3 below where the situation is discussed in the 
case of KdS geometry. Strong collimation of open equipotential surfaces near the axis of rotation suggests a tendency 
of backgrounds with A > to collimate streams of particles (jets) moving along the rotational axis. 

Equilibrium toroidal configurations of barotropic perfect fluid orbiting in the Kerr-de Sitter background were 
analyzed thoroughly in [12] where also the case of naked singularities was discussed. Here we present an overview of 
the basic results obtained in the approximation of test discs, and add an estimation for the central mass-energy density 
of an adiabatic fluid for which the obtained test-fluid results are valid. 



2. KERR-DE SITTER GEOMETRY 

The geometry of Kerr-de Sitter spacetimes is given by the line element 

d, 2 = -4h( dt - " sin2 S#) 2 + \ adt - (r 2 + « 2 ) d0l 2 + ^dr 2 + ^-dd 2 (4) 

I z p z I z p z 1 y ' J A r Ag 

where 

A r = r 2 -2Mr + a 2 -^Ar 2 (r 2 +a 2 ), A = 1 + ^Aa 2 cos 2 9 , I=l + ^Aa 2 , p 2 = r 2 + a 2 cos 2 0. (5) 

The spacetime is specified by three parameters: central mass (M), rotational parameter (a) corresponding to the 
specific angular momentum of the central object, and positive cosmological constant (A). It is convenient to introduce 
the dimensionless "cosmological parameter" 

y= l -AM 2 (6) 
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FIGURE 1. Division of KdS spacetimes according to the existence of stable equatorial circular orbits (SECOs). (I) Black-hole 
(BH) spacetimes with both direct and retrograde SECOs. (II) BH spacetimes with no retrograde SECOs. (Ill) BH spacetimes with 
no SECOs. (IV) Naked-singularity (NS) spacetimes with no SECOs. (V) NS spacetimes with both direct and retrograde SECOs. 
(VI) NS spacetimes with no retrograde SECOs. The dashed-dotted curve forms the boundary of NS region, where the retrograde 
plus-family SECOs exist; shaded is its subregion where such SECOs possess negative energy. 

and reformulate relations (4)-(5) into the completely dimensionless form by putting M = 1 hereafter. The spacetime 
is stationary, axially symmetric and asymptotically de Sitter. 

The spacetime horizons are determined by the condition A,. = 0. In general, three horizons (the inner and the outer 
of the black hole, r\- and r^ + , and the cosmological one, r c ) exist. A more detailed discussion on the existence of 
event horizons of KdS geometry including the extreme cases in which two or even three horizons coincide can be 
found in [7]. Note that cosmological horizon, behind which the geometry is dynamic, exists in all KdS spacetimes, 
and that there is a maximal value of the cosmological parameter allowing the existence of black holes: y cr it = 0.059. 

An analysis of the equatorial circular orbits of test particles in terms of the stability against radial perturbations 
enables to divide the parametric space (y, a 2 ) into six regions corresponding to the black-hole/naked-singularity 
spacetimes containing stable equatorial circular geodesies of a given family (Fig. 1). In the black-hole backgrounds, 
the plus-family corresponds to the stable corotating orbits while the minus-family to the stable counterrotating ones, 
where the direction of the orbits is related to the locally non-rotating frames (LNRFs) introduced by Bardeen et al. 
[13]; the orbits with locally measured azimuthal component of its 4-momentum > are called direct or corotating, 
while the ones with PW < are called retrograde or counterrotating. In the naked-singularity backgrounds, there is 
a sort of stable plus-family orbits which are counterrotating from the point of view of LNRFs, even with negative 
specific energy (the shaded region in Fig. 1). Note that also in black-hole backgrounds, contrary to the Kerr case 
(y = 0), the plus-family orbits can be both direct or retrograde, however, the retrograde ones, located near the upper 
limiting radius for the existence of the plus-family orbits (static radius or the retrograde photon circular orbit), are 
unstable. Moreover, for each family of orbits there is a maximal value of the cosmological parameter y allowing the 
stable circular orbits of a given family: y c ( ms+ ) = 0.069 and y c ( ms _) = 12/15 4 = 0.00024; the latter coincides with the 
maximal value of y allowing the stable circular orbits in SdS spacetimes [6]. A detailed discussion of the direction and 
stability of equatorial circular orbits can be found in [7]. 



Analytic theory of equilibrium configurations of rotating perfect fluid bodies was developed by Boyer [14] and than 
studied by many authors. The main result of the theory, known as "Boyer's condition", states that the boundary of any 
stationary, barotropic, perfect fluid body has to be an equipotential surface. Here we shall summarize its application 
to the relativistic test perfect fluid orbiting in a stationary and axisymmetric way in a stationary, axisymmetric 
background, as was introduced by Abramowicz and co-workers in the case of Schwarzschild and Kerr black holes 
[15, 16], and use it for an analysis in the KdS spacetimes [12]. 

In the standard Boyer-Lindquist coordinates the spacetime is described by the line element 



3. EQUIPOTENTIAL SURFACES 



ds 2 = g tt dt 2 + 2g t $dtd<j)+g H d<l) 2 +g rr dr 2 +g eg d0 2 



(7) 
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satisfying the properties of stationarity and axial symmetry, i.e., d t g^ v — d^g^y = 0. The stress-energy tensor of perfect 
fluid is given by 

T ti v = {e+p)U^U v +p8 v l , (8) 

where e and p are total energy density and pressure of the fluid, respectively. Further, we consider test perfect fluid 
with the 4-velocity = (V, , 0, 0); orbital motion of the fluid is characterized by vector fields of the angular 
velocity £l(r, 9) — /U' and the specific angular momentum £(r, 9) = —U^/U t , related by the metric coefficients 
of the background 

a = _8«+tg« _ (9) 



The relativistic Euler equation in the axially symmetric form reads 



where i = r, 9 and 

(7 7)2 _ 8%-gttgW 

[Ut) g„+2tg t i+Pg tt - (LL> 

For a barotropic fluid, i.e., for a body with an equation of state p = p(e), the surfaces of constant pressure are given 
by the equipotential surfaces of the potential W(r, 9) defined by the relations [16] 

dp :ln(I/ ( )m-ln(0i)+ f J^L. = Wia -W, (12) 



o £ + p v ' v ' Je^l-Sli 

where the subscript "in" refers to the inner edge of the disc. The explicit form of the potential, W = W(r, 9), is given 
by the relations (11) and (12), if one specifies the metric tensor of the background and the "rotational law", i.e., the 
function Q. = £1(1). The simplest but astrophysically very important is the case of uniform distribution of the specific 
angular momentum, 

£(r, 6) = const, (13) 

through the disc. It is well known that the tori with l(r, 9) = const are marginally stable [17] and capable to produce 
maximal luminosity at all [18]. Note that topological properties of the equipotential surfaces seem to be rather 
independent of the distribution of the specific angular momentum l(r, 9), see, e.g., [18, 19, 5]. In this special case 
the potential is given by very simple formula 

W(r,d)=lnU t {r,d). (14) 

The points where <9,W = correspond to free-particle (geodesic) motion due to the vanishing of the pressure-gradient 
forces there. Moreover, at the center of any perfect fluid torus the pressure attains the extreme value (maximum) and 
matter must follow a stable geodesic there. Thus, thick discs can exist only in the backgrounds allowing the motion 
along stable circular geodetical orbits. 

In the KdS backgrounds, the potential (14) takes the form 

Jp 2 A r A e sin 2 l' /2 

W(r, 6) - In | T - — + ^ - - ^ - _ - — | . (15) 

Performing the limit a — > 0, we get the corresponding form of the potential (14) in the SdS backgrounds [10]. All 
relevant properties of the equipotential surfaces are determined by the behaviour of the potential in the equatorial 
plane (9 = 7t/2). Inspecting the reality conditions of W(r, 9 =7t/2) in the stationary parts of the background (A r > 0; 
in the black-hole backgrounds, we restrict our attention to the stationary region between the (outer) black-hole and 
cosmological horizon only), we arrive at the condition for the occurrence of matter with a given distribution of l(r, 9), 

e P h-<£<e P h+, (16) 

where 

r 2 

£ph±(r,a,y)=a+ (17) 
a ± v A- 



Basic properties of toroidal structures in Kerr-de Sitter backgrounds 



February 5, 2008 



4 



3 




(log r) sin 9 (log r) sin (log r) sin 9 

(a) (b) (c) 

FIGURE 2. Typical behaviour of equipotential surfaces (meridional sections) in SdS and KdS black-hole spacetimes. Light gray 
region contains closed equipotential surfaces. The last closed surface is self-crossing in the cusp(s). Possible toroidal configurations 
correspond to: (a) accretion discs, (b) marginally bound accretion discs and (c) excretion discs. (The figures are drawn for the KdS 
spacetime with y = 1CT 6 , a 2 = 0.6 using the standard Boyer-Linquist coordinates.) 



corresponds to the effective potential of the photon geodesic motion; see [22] for an alternative definition. Local 
extrema of the function W(r, = n/2) lie at those radii where the specific angular momentum coincides with the 
specific angular momentum of test particles moving on the geodetical (Keplerian) circular orbits, i.e., where 

(r 2 + a 2 )(l-yry/ 2 T ar l ^ + r(r 2 + a 2 )y} 

lX±Y,a,y)-zt r 3/2[ 1 _( r 2 + fl 2 )y ]_ 2r l/2 ±a ( 1 _ };r 3)l/2 ' ^> 

Those extrema are the only local extrema of the function W(r, ). 

Toroidal configurations arise for such a distribution of l(r, 8) in the disc which intersects the Keplerian distribution 
£ K± (r) in the part(s) corresponding to stable circular orbits. In black-hole backgrounds, stationary toroidal configu- 
rations exist for I e (4ns(i): i ms ( )), where 4ns(i) (Ans(o)) corresponds to the Keplerian specific angular momentum on 
the inner (outer) marginally stable orbit. The same is true also in most of the naked-singularity backgrounds, however, 
exceptions exist concerning the plus-family discs in naked-singularity backgrounds with the rotational parameter low 
enough to admit counterrotating stable plus-family circular geodesies, see Fig. 1. In fact, there are naked-singularity 
backgrounds, in which the stationary tori exist for any uniform distribution of £(r, 9) in the disc [12]. 

In both black-hole and naked-singularity backgrounds we can distinguish three kinds of discs (Figure 2): 

accretion discs: Toroidal equipotential surfaces are bounded by the marginally closed critical equipotential surface 
self-crossing in the inner cusp and enabling outflow of matter from the disc into the BH/NS. Another critical 
surface self-crossing in the outer cusp is open. Matter filling the region between the critical surfaces cannot 
remain in hydrostatic equilibrium and contributes to the accretion flow along the inner cusp and a throat formed 
by open surfaces. Moreover, if the potential levels corresponding to the critical surfaces are comparable, i.e., 
Wcrit(i) < Wcrit(o)> huge overfilling of the critical surface with the inner cusp causing the accretion could be 
combined with the so-called excretion, i.e., outflow through the outer cusp (after overfilling of the critical surface 
with the outer cusp), having a capability to regulate the accretion. 

marginally bound accretion discs: Such configurations exist only for the uniform distribution of the specific angular 
momentum in the disc £(r, 0) = £ m b, where t m \, corresponds to the Keplerian specific angular momentum on the 
marginally bound circular orbit. Toroidal equipotential surfaces are bounded by the marginally closed critical 
equipotential surface self-crossing in both the cusps. Any overfilling of the critical surface causes the accretion 
inflow through the inner cusp as well as the excretion outflow through the outer cusp. 

excretion discs: Toroidal equipotential surfaces are bounded by the marginally closed critical equipotential surface 
self-crossing in the outer cusp and enabling outflow of matter from the disc into the outer space by a violation 
of hydrostatic equilibrium. The equipotential surface with the inner cusp, if such a surface exists, is open 
(cylindrical) and separated from the critical surface with the outer cusp by other cylindrical surfaces which, 
in fact, disable accretion into the black hole. 
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FIGURE 3. Examples of "exotic" toroidal configurations (meridional sections through the equipotential surfaces) in KdS naked- 
singularity spacetimes. Shaded regions contain closed equipotential surfaces. The last closed surface is self-crossing in the cusp(s). 
(a) Counterrotating negative-energy accretion disc, (b) Two isolated discs - the inner one is a counterrotating accretion disc, the 
outer one is a corotating excretion disc, (c) Two isolated discs - the inner one is a counterrotating accretion disc, the outer one is a 
marginally bound counterrotating accretion disc. (The figures are drawn for NS spacetimes with y = 10~ 6 , a 2 = 1 .05 and the values 
of £ being equal to 8, 3.5 and —4.72, respectively. The Kerr-Schild coordinates, covering the whole range of the disc including the 
ring singularity, and their appropriate scaling are used.) 

In KdS naked-singularity backgrounds allowing counterrotating stable plus-family circular orbits (Fig. 1), some 
more exotic configurations exist: 

a. counterrotating negative-energy discs 

Specific energy of the fluid elements in the center and the cusps (where the fluid follows the geodesic motion) is 
negative and we can expect that every fluid element in the disc has energy E < 0. Moreover, no open equipotential 
surfaces going out from the singularity are connected with such configurations. In the case of accretion discs, the 
equipotential surface with the outer cusp need not exist (Fig. 3a). 

b. configurations with two isolated discs corresponding to the same value of £ 

The inner disc is always counterrotating accretion disc (with matter in the states with E < or E > 0), the outer 
disc can be, in dependence on the value of £, the corotating or counterrotating excretion disc, as well as the 
counterrotating accretion disc (Figs. 3b, c). 

Due to the existence of non-zero pressure-gradients in the fluid, the inner edge of accretion discs (corresponding 
to the inner cusp of equipotential surfaces) is shifted under the inner marginally stable circular orbit up to the inner 
marginally bound circular orbit, r mb ^ < r m < r ms ^y Similarly, the outer edge of excretion discs (corresponding to 
the outer cusp of equipotential surfaces) is located between the outer marginally stable and outer marginally bound 
circular orbit, r ms ( ) < r out < r mb / y Marginally bound accretion discs have, thus, naturally determined both edges by 
the location of the cusps of the only critical surface, r m r mh ^, r out r mb ( ), and correspond to maximally extended 
discs. Moreover, potential difference between the boundary (determined by the marginally closed critical surface) 
and the center of the torus, AW = W CI \t — Winter, takes the largest values for plus-family marginally bound accretion 
discs. In black-hole backgrounds, the maximal value corresponds to the disc corotating the extreme Kerr black hole 
(y = 0), AW i=a 0.549 [16], and with the cosmological parameter y growing up to y cr it = 0.059 tends to zero. In naked- 
singularity backgrounds, the potential difference grows unlimitedly, AW — ► °°, for the plus-family discs orbiting a 
naked singularity approaching the extreme-hole state, independently of the cosmological parameter y < y cr j t . 



4. MASS ESTIMATES OF ADIABATIC TORI 



The total mass of the disc, m, is given by the Tolman's formula [15] 




(19) 
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TABLE 1. (Upper part) Mass parameter, the static radius and radius of the outer 
marginally stable orbit in extreme KdS black-hole spacetimes with the current value of 
the cosmological constant, Ao ~ 1.3 x 10~ 56 cm~ 2 . (Lower part) Central mass-density 
of an adiabatic fluid for which the total mass of the disc is comparable with the mass 
of the SdS black hole for two values of adiabatic index y = 5/3 and 7=7/5. 
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where g — det(g^ v ) and Ty is the stress-energy tensor of perfect fluid given by the relation (8). Assuming an adiabatic 
equation of state for a barotropic perfect fluid [24] 

P = K P y, y=i + - (20) 

n 

where 7 is an adiabatic index, and a non-relativistic limit with p<^e~ p, where p is the rest-mass density, the relation 
(19) can be written in an approximative form [15] 

m = 2np c I 

./disc 

in which the functions £l(r, 9) and W(r, 9) are given by the relations (9) and (15). W c and p c correspond to the potential 
value and the rest-mass density in the centre of the disc, respectively. 

Comparing the total mass of the disc m with the mass of the black hole M, i.e. m w M, we can get the maximal value 
of the central mass-density for which the approximation of test fluid is valid in the sense m <C M. The first results of 
numerical integration have been obtained for marginally bound accretion discs orbiting in the SdS background and for 
two principal values of adiabatic index 7=5/3 and 7 = 7/5. The results are presented in Table 1 . 



l+£Cl(r,d) 

i-eci(r,e) 



W; n -W c 



(r 2 + a 2 cos 2 0)sin0drd0 



(21) 



5. CONCLUSIONS 

We conclude that there are two features in the structure of equipotential surfaces near the static radius of a given space- 
time connected with a cosmic repulsion; both have already been known from the previous analysis of equipotential 
surfaces in SdS backgrounds [10]. The first one is the outer cusp which enables the outflow of matter from the new type 
of stationary tori — the excretion discs. However, the outer cusp plays an important role also for the accretion discs, as 
puts an upper limit on their extension. Recall that no such an outer edge is naturally defined in the case of accretion 
discs orbiting a single black hole (naked singularity) in asymptotically flat Schwarzschild or Kerr backgrounds. The 
second feature connected with the cosmic repulsion consists in strong collimation of open equipotential surfaces near 
the axis of rotation, being evident nearby and behind the static radius, when compared with the Kerr case, see [12] for 
an illustrative figure, suggesting a certain role of A > in the collimation of jets far away from the maternal accretion 
or excretion disc. Rotation of the background influences the shape of tori: the corotating discs are thicker and more 
extended than the discs in corresponding SdS background, generating a narrower funnel where the jets are most prob- 
ably created [25]. On the other hand, the same is true for the discs in SdS backgrounds when these are compared with 
counterrotating discs in corresponding KdS background. 

Finally, we shall give an idea on scales, at which the discussed effects take place, by expressing basic characteristics 
of the tori in astrophysical units for the current value of the cosmological constant, A = Ao. The results are presented 
in Table 1. As the outer edge of tori is located between the outer marginally stable orbit and the static radius, 
r ms(o) < r out < r s , the repulsive cosmological constant puts a limit on maximal extension of disc-like structures in a 
given background. Remarkably for supermassive black holes (1O 6 M -1O 1O M ), the dimensions of test tori are roughly 
comparable with the dimensions of galaxies [26]. The upper limits on the central mass-density of the adiabatic tori 
orbiting supermassive SdS black holes (10 8 -10 9 M Q ), for which the test-fluid approximation is valid, are in the case 
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of adiabatic index y = 7/5 one or two orders higher than the typical density of molecular clouds 10~ 18 kgrrT 3 [26]. 
Of course, to get a more realistic picture, influence of the cosmic repulsion on self-gravitating tori has to be studied. 
On the other hand, since the jets escaping from AGN can many times exceed the dimension of the galaxy, repulsive 
cosmological constant could play an important role in the collimation of jets far away from the "seed" galaxy. 
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